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Abstract — The paper proposes a framework for modeling and 
analysis of the dynamics of supply, demand, and clearing prices 
in power system with real-time retail pricing and information 
asymmetry. Real-time retail pricing is characterized by passing 
on the real-time wholesale electricity prices to the end consumers, 
and is shown to create a closed-loop feedback system between the 
physical layer and the market layer of the power system. In the 
absence of a carefully designed control law, such direct feedback 
between the two layers could increase volatility and lower the 
system's robustness to uncertainty in demand and generation. A 
new notion of generalized price-elasticity is introduced, and it is 
shown that price volatility can be characterized in terms of the 
system's maximal relative price elasticity, defined as the maximal 
ratio of the generalized price-elasticity of consumers to that of the 
producers. As this ratio increases, the system becomes more volatile, 
and eventually, unstable. As new demand response technologies 
and distributed storage increase the price-elasticity of demand, the 
architecture under examination is likely to lead to increased volatility 
and possibly instability. This highlights the need for assessing 
architecture systematically and in advance, in order to optimally 
strike the trade-offs between volatility, economic efficiency, and 
system reliability. 

Index Terms — Real-Time Pricing, Volatility, Lyapunov Analysis. 
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I. Introduction 

THE increasing demand for energy along with growing en- 
vironmental concerns have led to a national agenda for 
engineering a modem power grid with the capacity to integrate 
the renewable energy resources at large scale. In this paradigm 
shift, demand response and dynamic pricing are often promoted as 
means of mitigating the uncertainties of the renewable generation 
and improving the system's efficiency with respect to economic 
and environmental metrics. The idea is to allow the consumers 
to react-in their own monetary or environmental interest-to 
the wholesale market conditions, possibly the real-time prices. 
However, this real-time or near real-time coupling between supply 
and demand creates new challenges for power system operation. 
The source of a most significant challenge is the information 
asymmetry between the consumers and the system operators. 
Indeed, real-time pricing under information asymmetry induces 
additional uncertainties due to the uncertainty in consumer behav- 
ior, preferences, private valuation for electricity, and consequently, 
unpredictable reactions to real-time prices. 

The existing body of literature on dynamic pricing in commu- 
nication or transportation networks is extensive. See for instance 
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ifTTl . (6), lfl4l and the references therein. However, the specific 
characteristics of power systems which can be attributed to 
uncertainty in consumer behavior, the close coupling and real- 
time interaction of economics and physics, and the reliability and 
operational requirements that supply must match demand at all 
times raise very unique challenges that need to be addressed. 

Various forms of dynamic retail pricing of electricity have been 
advocated in economic and engineering texts. In J5j, Borenstein 
et. al. study both the theoretical and the practical implications of 
different dynamic pricing schemes such as Critical Peak Pricing, 
Time-of-Use Pricing, and Real-Time Pricing. They argue in favor 
of real-time pricing, characterized by passing on a price, that 
best reflects the wholesale market prices, to the end consumers. 
They conclude that real-time pricing delivers the most benefits 
in the sense of reducing the peak demand and flattening the load 
curve. In ifTOl . Hogan identifies dynamic pricing, particularly real- 
time pricing as a priority for implementation of demand response 
in organized wholesale energy markets. Similar conclusions are 
reached in a study conducted by Energy Futures Australia fl8] . 

The appeal of dynamic retail pricing is not limited to theoretical 
research and academic studies, and real-world implementations 
are emerging at a rapid pace. For instance, California's state's 
Public Utility Commission has enacted a series of new energy 
regulations which set a deadline of 2011 for the state utilities to 
propose a new dynamic pricing rate structure, specifically defined 
as an electric rate structure that reflects the actual wholesale 
market conditions, such as critical peak pricing or real-time 
pricing ll24l . In this paper, we show that directly linking the 
consumer prices to the wholesale market prices creates a close- 
loop feedback system with the Locational Marginal Prices as 
the state variables. We observe that such feedback mechanisms 
may increase volatility and decrease the market's robustness to 
uncertainty in demand and generation. We introduce a notion of 
generalized price-elasticity, and show that price volatility can be 
upperbounded by a function of the system's Maximal Relative 
Price-Elasticity (MRPE), defined as the maximal ratio of the 
generalized price-elasticity of consumers to the generalized price- 
elasticity of producers. As this ratio increases, the system may 
become more volatile, eventually becoming unstable as the MRPE 
exceeds one. 

While the system can be stabilized and volatility can be 
reduced in many different ways, e.g., via static or dynamic 
controllers regulating the interaction of wholesale markets and 
retail consumers, different pricing mechanisms pose different 
consequences on competing factors of interest such as volatility, 
operational reliability, economic efficiency, and environmental 
efficiency. The intended message is that the design of a real- 
time pricing mechanism must take system stability issues into 
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consideration, and that successful design and implementation of 
such a mechanism entails careful analysis of consumer behavior 
in response to price signals, and the trade-offs between volatility, 
reliability, and economic or environmental efficiency. 

Prior research relevant to stability of power markets has 
appeared in several papers by Alvarado [fl], lEI on dynamic 
modeling and stability, Watts and Alvarado [4| on the influence 
of future markets on price stability, and Nutaro and Protopopescu 
|[T5l on the impact of market clearing time and price signal 
delay on power market stability. The model adopted in this 
paper differs from those of 0], E], 0, and lfT31 in that we 
analyze the global properties of the full non-linear model as 
opposed to the first-order linear differential equations examined 
in these papers. In addition, the price updates in our paper 
occur at discrete time intervals and are an outcome of marginal 
cost pricing in the wholesale market by an Independent System 
Operator (ISO), which is consistent with the current practice 
in deregulated electricity markets. Furthermore, beyond stability, 
we are interested in providing a characterization of the impacts 
of uncertainty in consumer behavior on price volatility and the 
system's robustness to uncertainties. 

The organization of this paper is as follows. In Section UT1 we 
present some preliminary concepts and definitions. In Section Hill 
we present a mathematical model for the dynamic evolution of 
supply, demand, and clearing prices under real-time pricing. Sec- 
tion [IV] contains the main theoretical contributions of this paper: 
characterizing volatility in terms of the market's maximal relative 
elasticity and uncertainty in consumer behavior. In Section [V] 
we qualitatively discuss our results, compare with some of the 
results in the literature, and point to some important questions 
regarding the trade-offs arising due to uncertainty in generation 
and quantifying the value of information. Numerical simulations 
are presented in Section [VT] Finally, we offer some closing 
remarks and further directions for future research in Section [VII] 

II. Preliminaries 

A. Notation 

The set of positive real numbers (integers) is denoted by R + 
(Z + ), and nonnegative real numbers (integers) by R + (Z+). The 
class of real-valued functions with a continuous n-th derivative on 
IcMis denoted by C n X. For a vector v eM. l ,Vk denotes the fc- 

def 

th element of v, and ||u|L denotes the standard p-norm: ||u|L = 
/ , \ i/p 

[J2i=i \ v i\ ) ■ Also, we will use ||y|| to denote any p-norm 
when there is no ambiguity. The space of M 1 -valued functions 
h : Z i-> R ( of finite p-norm 

\\m= E \\Ht)\%= E EiM^r 

t— — oo t— — oo i—l 

is denoted by i p (Z) or simply l v when there is no ambiguity. 
For a differentiable function / : R" H> R m , we use / to denote 
the Jacobian matrix of /. When / is a scalar function of a single 
variable, / simply denotes the derivative of / with respect to 
its argument: / (x) = df(x)/dx. Since throughout the paper 



time is a discrete variable, this notation would not be confused 
with derivative with respect to time. Finally, for a measurable set 
X C R, /xl (X) is the Lebesgue measure of X. 

B. Basic Definitions 

1) Volatility: 

Definition 1: Scaled Incremental Mean Volatility (IMV): 

Given a sig nal h : Z i-> R { , and a function p : M. 1 h-» R m , the 
p-scaled incremental mean volatility measure of h (■) is defined 
as 

1 T 

V p (h) = lim -E IIP (h (t + 1)) - p (h (t))\\ (1) 

t=0 

where, to simplify the notation, the dependence of the measure 
on the norm used in (HJ is dropped from the notation V p (h). 

To quantify volatility for fast-decaying signals with zero or 
small scaled IMV, e.g., state variables of a strictly stable au- 
tonomous system, we will use the notion of scaled aggregate 
volatility, defined as follows. 

Definition 2: Scaled Incremental Aggregate Volatility 
(IAV): Given a signal /i:Z^l ! , and a function p : R z i-> R"\ 
the p-scaled incremental aggregate volatility measure of h (•) is 
defined as 

oo 

V p (^) = ^ + 1)) -p(/i(i))|| . (2) 

t=0 

In particular, we will be interested in the log-scaled incremental 
volatility as a metric for quantifying volatility of price, supply, or 
demand in electricity markets. 

Remark 1: The notions of incremental volatility presented in 
Definitions Q] and [2] accentuate the fast time scale, i.e., high fre- 
quency characteristics of the signal of interest. Roughly speaking, 
the scaled IMV or IAV are measures of the mean deviations of 
the signal from its moving average. In contrast, sample variance 
or CV (coefficient of variation, i.e., the ratio of standard deviation 
to mean) provide a measure of the mean deviations of the 
signal from its average, without necessarily emphasizing the high- 
frequency characteristics. A slowly-varying signal with a large 
dynamic range may have a large sample variance or CV, but a 
small IMV, and thus will be considered less volatile than a fast- 
varying signal with a large scaled IMV. Since we are interested 
in studying the fast dynamics of spot prices in electricity markets 
and the associated stability/reliability threats, the scaled IMV and 
IAV as defined above are more appropriate measures of volatility 
than variance or CV. 

2) Stability: The notion of stability used in this paper is the 
standard notion of asymptotic stability. Consider the system 

x{t+l)=ij){x{t)) (3) 

where ip (■) is an arbitrary map from a domain X C R™ to R™. 
The equilibrium i e X of (01 is stable in the sense of Lyapunov 
if all trajectories that start sufficiently close to x remain arbitrarily 
close to it, i.e., for every e > 0, there exists S > such that 

\\x (0) - x\\ < S ^ \\x (t) - x\\ < e, Vt>0 

The equilibrium is globally asymptotically stable if it is Lyapunov 
stable and for all x (0) G X : lim^oo x (t) — x. 
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C. Market Structure 

We begin with developing an electricity market model with 
three participants: 1. The suppliers, 2. The consumers, and 3. 
An Independent System Operator (ISO). The suppliers and the 
consumers are price-taking, profit-maximizing agents. The ISO 
is an independent, profit-neutral player in charge of clearing 
the market, that is, matching supply and demand subject to the 
network constraints with the objective of maximizing the social 
welfare. Below, we describe the characteristics of the players in 
more detail. 

1) The Consumers and the Producers: Let D = {1, . . . ,n g } 
and S = {1, . . . , n s } denote the sets of consumers and producers 
respectively. Each consumer j G D is associated with a value 
function Vj : M.+ n- E, which can be thought of as the monetary 
value that consumer j derives from consuming x units of the 
resource, electricity in this case. Similarly, each producer i G S, 
is associated with a function a : R+ i- > K+ representing the 
monetary cost of production of the resource. 

Assumption I: For all i £ S, the cost functions c,; (•) are in 
C 2 (0, oo), strictly increasing, and strictly convex. For all j G D, 
the value functions Vj (•) are in C 2 (0, oo), strictly increasing, and 
strictly concave. 



Let <L 



j G D, and 



i £ S be. 



demand and supply functions mapping price to consumption and 
production, respectively. In the framework adopted in this paper, 
the producers and consumers are price-taking, utility-maximizing 
agents. Therefore, letting A be the price per unit of electricity, we 
have 

dj (A) = arg max Vj (x) — Ax, j G -D, (4) 



max{0, {x | ij (x) = A}} 



and 



Si (A) — arg max Xx — Ci (x) , i G S. (5) 

= max{0, {x | Ci (x) = A}} 

For the sake of convenience in notation and in order to avoid un- 
necessary technicalities, unless stated otherwise, we will assume 
in the remainder of this paper that dj (A) = vj (A) is the demand 
function, and Si (A) = (A) is the supply function. This can 
be mathematically justified by assuming that v (0) = oo, and 
c (0) = 0, or that A G [c (0) , ii (0)] . 

Definition 3: The social welfare S is the aggregate benefit of 
the producers and the consumers: 

jeD zes 

If A; = Aj = A, Vi, j, we say that A is a uniform market clearing 
price, and in this case, we have: 

j£D ieS 



a) Heterogeneous Consumers with Uncertain Value Func- 
tions: We will consider two models of heterogenous consumers 
with time-varying value functions to represent the uncertainty in 
consumer behavior. 

- Multiplicative Perturbation Model: The uncertainty in 
consumer's value function is modeled as 



Vj (x, t) = aj (t) v 



a 3 (i) J ' 



J G D, 



(6) 



where v a : R + H» M is a nominal value function and aj : Z + M- 
R + is an exogenous signal or disturbance. Given a price A (t) > 0, 
under the multiplicative perturbation model (|6) we have 



dj(X,t) =a.j(t) Vq 1 (A (*)) 



(7) 



Thus, the same price A may induce different consumptions at 
different times, depending on the type and composition of the 
load. 

- Additive Perturbation Model: The uncertainty in con- 
sumer's value function is modeled as 



ij (x, t) =v {x- Uj (t)) , jeD, 



(8) 



where u, 



f h-> R+ is exogenous. Thus, given a price A (t) > 0, 
under the additive perturbation model (|8j, the demand function 
is 



d j (X,t)=u j (t) + v- 1 (X(t)) 



(9) 



- Aggregation of Several Consumers: The aggregate re- 
sponse of several consumers (or producers) to a price signal 
may be modeled as the response of a single representative agent, 
although explicit formula for the utility of the representative agent 
may sometimes be too complicated to find flip , ID- F° r the case 
of N identical consumers with value functions Vj = v , j G D, 
it can be verified that the aggregate demand is equivalent to the 
demand of a representative consumer with value function ITT81 : 



v(x) = Nv (^ 



(10) 



Suppose now, that the consumer behavior can be modeled via 
©-(01. Let 



N 



« («) = 2^ J=1 a i (*) . 



and suppose that there exists a nominal value ao, such that 

a (t) =a + Aa (t) = a (1 + S (*)) 

where 6 (t) — Ace(t) /ao satisfies \S (t)\ < 1. Define v (x) — 
ctQV (x/ao) ■ It can be then verified that the aggregate demand 
can be modeled as the response of a representative agent with 
value function 

v (x, t) = a (t) v 



a(t) 
(a + Act (t)) v 

(1 + S(t))v 



a + Aa (t) 



l + 5(t) 



(11) 
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The aggregate response is then given by 

d(A(t),t) = (l + Sit))^ 1 (\(t)) 



(12) 



Similarly, under the additive perturbation model the aggregate 
behavior can be represented by 



v (x, t) 
d(X (t),t) 



v [x - 
«(*) 



«(*)) 



(13) 
(14) 



where v (•) is given by < TT~Qb and tt (i) = Uj (t) . The interpre- 
tation of ( foi l and (fl4l is that at any given time t, the demand 
comprises of an inelastic component u (t) which is exogenous, 
and an elastic component v^ 1 (X(t)). Another interpretation is 
that v -1 (A (t)) represents the demand of those consumers who 
are subject to real-time pricing, and u (t) represents the demand 
of the non-participating consumers. 

2) The Independent System Operator (ISO): The ISO is a non- 
for-profit entity whose primary function is to optimally match sup- 
ply and demand subject to network and operational constraints. 
The constraints include power flow constraints (Kirchhoff 's laws), 
transmission line constraints, generator capacity constraints, local 
and system-wide reserve capacity requirements and possibly some 
other constraints specific to the ISO ll23l . lETI . Il22l . For real- 
time market operation, the constraints are linearized near the 
steady-state operating point and the ISO optimization problem is 
reduced to a convex-typically linear-optimization often referred 
to as the Economic Dispatch Problem (EDP), or the Optimal 
Power Flow Problem. A set of Locational Marginal Prices (LMP) 
emerge as the dual variables corresponding to the nodal power 
balance constraints. These prices vary from location to location 
as they represent the marginal cost of supplying electricity to a 
particular location. We refer the interested reader to [E01 . Il22l . 
Il2p . and ifTTl for more details. However, we emphasize that the 
spatial variation in the LMPs is a consequence of congestion 
in the transmission lines. When there is sufficient transmission 
capacity in the network, a uniform price will materialize for the 
entire system. With this observation in sight, and in order to 
develop tractable models that effectively highlight the impacts of 
the behavior of producers and consumers-quantified through their 
cost and value functions-on system stability and price volatility, 
we will make the following assumptions: 

1) Resistive losses are negligible. 

2) The line capacities are high enough, so, congestion will not 
occur. 

3) There are no generator capacity constraints. 

4) The system always has sufficient reserve capacity and the 
marginal cost of reserve is the same as the marginal cost 
of generation. 

Under the first two assumptions, the network parameters be- 
come irrelevant in the supply-demand optimal matching problem. 
The third and fourth assumptions are made in the interest of 
keeping the development in this paper focused. They could, 
otherwise, be relaxed at the expense of a somewhat more involved 
technical analysis. The last assumption also implies that we do not 
differentiate between actual generation and reserve. A thorough 
investigation of the effects of network constraints and reserve 



capacity markets, whether they are stabilizing or destabilizing, 
does not fall within the scope of this paper. The interested readers 
may consult Q, IfTTl for an analysis of dynamic pricing in 
electricity networks with transmission line and generator capacity 
constraints. 

Under the above assumptions, the following problem charac- 
terizes the ISO's optimization problem: 



max 2J v j (dj ) - 2jcj(sj 

j£D ies 



s.t. 



J2 d J = & 

jeD ies 



(15) 



The following lemma which is adopted from [12], provides the 
justification for defining the LMPs as the Lagrangian multipliers 
corresponding to the balance constraint(s). 

Lemma 1: Let d* — [d\, ■ • • , d* ] , and s* — [s*, ■ ■ • , s*J 
where d*, j G D and s* , i G S, solve (Q3). There exists 
a price A* G (0, oo) , such that d* and s* solve © and ©. 
Furthermore, A* is the Lagrangian multiplier corresponding to 
the balance constraint in (15[ . 

Proof: The proof is based on Lagrangian duality and is 
omitted for brevity. The proof in |fi2) would be applicable here 
with some minor adjustments. ■ 

The implication of Lemma Q] is that by defining the market 
price to be the vector of Lagrangian multiplier corresponding to 
the balance constraints, the system operator creates a competitive 
environment in which, the collective selfish behavior of the 
participants results in a system-wide optimal condition. 

a) Real-Time System Operation and Market Clearing: 
Consider the case of real-time market operation and assume that 
price-sensitive retail consumers do not bid in the real-time market. 
In other words, they do not provide their value functions to the 
system operator (or any intermediary entity in charge of real- 
time pricing). Though, they may adjust their consumption in 
response to a price signal, which is assumed in this paper, to 
be the wholesale market clearing price. In this case, the demand 
is assumed to be inelastic over each short pricing interval, and 
supply is matched to demand. Therefore, (fT5t reduces to meeting 
the fixed demand at minimum cost: 



mm y^Cjjsj 

ies 

s.t. J2 s * = 5^ 

tes jeD 



(16) 



where dj is the predicted demand of consumer j for the next time 
period. We assume that the system operator solves (fTSl and sets 
the price to the marginal cost of production at the minimum cost 
solution. The discrepancy between scheduled generation (which is 
equal to the predicted demand) and actual demand is compensated 
through reserves with the same marginal costs. Thus, we will not 
include reserve parameters and equations explicitly in the model. 
More details regarding a dynamic extensions of this model are 
presented in the next Subsection. 
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s(t + 1) = d(t) 

X(t+l)=^c(x) 
ax . 



d(t) = argmax v(x) — X(t)x 



Fig. 1. Exante Priced Supply/Demand Feedback 



III. Dynamic Models of Supply-Demand under 
Information Asymmetry 

In this section, we develop dynamical system models for the 
interaction of wholesale supply and retail demand in electricity 
markets with information asymmetry. In this context, "asymmetry 
of information" refers to the architecture of the information layer 
of the market, in which, the market operator has full information 
about the cost of supplying the resource (e.g., through the offers 
of the producers), but has no information about valuation of the 
resource by the demand side. 

The real-time market is cleared at discrete time intervals and 
the prices are calculated and announced for each interval The 
practice of defining the clearing price corresponding to each 
pricing interval based on the predicted demand at the beginning 
of the interval is called exante pricing. As opposed to this, ex- 
post pricing refers to the practice of defining the clearing price 
for each pricing interval based on the materialized consumption 
at the end of the interval. In ex-post pricing the demand is subject 
to some price uncertainty as the actual price will be revealed after 
consumption has materialized. In exante pricing without ex-post 
adjustments, the entity in charge of real-time pricing faces the 
price uncertainty!!, as it will have to reimburse the generators 
based on the actual marginal cost of production, while it can 
charge the demand only based on the exante price. We will present 
dynamic market models for both pricing schemes. These models 
are consistent with the current practice of marginal cost pricing 
in wholesale electricity markets, with the additional feature that 
the retail consumers adjust their usage based on the real-time 
wholesale market price. 

A. Price Dynamics under Exante Pricing 

Let A (t) denote the exante price corresponding to the consump- 
tion of one unit of electricity in the time interval [t, t + 1] . Let 
d{t) — J2jeD^i (*) b e the actual aggregate consumption during 

'in most regions of the United States, such as New England, California, or 
PJM, the real-time market is operated in five-minute intervals. 

2 ln this paper we combine the role of the ISO and that of an entity in charge 
of real-time pricing in the retail market. Whether in practice this will be the case 
or not, has no influence on the intended message and the results that we deliver. 



this interval: 

d (t) = = 5>7 X (A (t)). (17) 

jeD jeD 

Since Vj (•) is known only to consumer j, at time t, only an 
estimate of d (t) is available to the ISO, based on which, the price 
A (t) is calculated. The price A (t) is therefore, the marginal cost 
of predicted supply to meet the predicted demand for the time in- 
terval [t, We assume that the ISO's predicted demand/supply 
for each time interval ahead is based on the actual demand at the 
previous intervals: s t +i = d (i + 1) = <j> (d(t) , ■ ■ ■ ,d(t — T)) , 
T € Z. The following equations describe the dynamics of the 
market: 

cr 1 ( A (* + 1)) = * (* + 1) = d{t + 1) (is) 

ies 

d(i+l) = 0(d(t),"- ,d(t-T)) (19) 
vj 1 (A (t-k)) = d(t-k), Vfc < T (20) 

jeD 

where d20]l follows from (H3, and A (i+1) in (H) is the 
Lagrangian multiplier associated with the balance constraint 
in optimization problem ( TToT l solved at time t + 1, i.e., with 
Y,^ D dj = d{t + l). 

The prediction step ( fl9l l may be carried through by resorting 
to linear auto-regressive models, in which case, we will have: 

T 

<f>(d{t),--- ,d(t-T)) = J2<*kd(t-k), a k eR. (21) 

fc=0 

When (f>(-) is of the form 1211 . equations (Ti8T> — (f20b result in: 

T 

J2 c- 1 (a (t +i)) = Y, a *Y, v 1 ( A - fc )) ( 22 ) 

ies k=a jeD 

Some of the popular forecasting models currently in use by the 
system operators are variations of the persistence model which 
corresponds to the special case where the predicted demand for 
the next time step is assumed to be equal to the demand at the 
previous time step, i.e., <j){d[t) , • • • , d (t — T)) = d (t). In this 
case, equations (fT~8b — (f20t result in: 

^c,- 1 (A(t + l)) = ^«- 1 (A(<)). (23) 

ies jeD 

If all the producers can be aggregated into one representative 
producer agent with a convex cost function c(-), and all the 
consumers can be aggregated into one representative consumer 
agent with a concave value function v (•) , then d22b and (123b 
reduce, respectively, to : 

\(t + l) = d(^2 T k=Q a k v~ 1 (\(t-k)^j (24) 

and 

X(t + l) = c:(i- 1 (X(t))). (25) 

More details on the construction of the representative agent 
mode can be found in fl8l . 
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s(t) = d(t) 




A(t)=|e(.) 


.(*} 



A(t + l)=A(t) 

d(t + l)=argmaxti(xl— X(t + lis 

Fig. 2. Ex-post Priced Supply/Demand Feedback 



B. Price Dynamics under Ex-post Pricing 

Under ex-post pricing, the price charged for consumption of 
one unit of electricity during the interval [t, t + 1] is declared at 
time t+1, when the total consumption has materialized. In order 
to decide on the amount to consume during each time interval 
ahead, a prediction of the ex-post price is needed. We assume 
that Xj (t + 1) , consumer j's predicted price, is a function of the 
ex-post prices of the previous intervals. Therefore, 

X j (t + l) = (f> j (X(t),--- ,A(t-T)) (26) 

rf ( i + 1 ) = E ,) 7 1 ( AJ ( < + 1 )) ( 27 > 

jGD 

J2K 1 (*(t+l)) = d(t+l) (28) 

ies 

By combining (f26b — (f28b we obtain: 

(A (t + !)) = £ vj 1 (<f>j (A (i) , • • • , A (t - T))) (29) 

It is observed that when the consumers use the per- 
sistence model for predicting future prices, i.e., when 
(j>j (A (t) , ■ ■ ■ , A (t — T)) = A (i) , Vj, then the price dynamics 
d29l becomes identical to the case with exante pricing (1231 , with 
the difference that the price uncertainty and the associated risks 
are bore by the consumer. In general, however, the price dynamics 
would depend on how each individual consumer predicts the 
ex-post price. This additional layer of dependency on consumer 
behavior suggests that more complicated market outcomes with 
multiple, possibly inefficient equilibria could materialize in ex- 
post-priced retail markets. 

Remark 2: Equation (f28b assumes that the generators were 
dispatched optimally, which is ideal but unlikely in practice. 
In this paper, we do not model the intricacies arising from the 
discrepancy between exante dispatch (which is the actual dispatch 
schedule based on prediction, and hence, not necessarily optimal) 
and ex-post dispatch (which characterizes how the generators 
should have been ideally dispatched). Although the settlement 
of these discrepancies is important in practice, such details can 
be safely ignored without affecting the core of our framework. 



Remark 3: It is also possible to consider dynamic models aris- 
ing from exante pricing complemented with ex-post adjustments, 
see for instance |fl9l . 

C. Demand Dynamics under Exante or Ex-post Pricing 

We could alternatively write dynamical system equations for 
the evolution of demand. Under exante pricing we will have: 

vj {dj (t + 1)) = c, (m (t)) Vz G S, j e D (30) 

E s > (*) - * (E ieD d i (*)>-' E ieD d i (« - T )) < 

ies 

(31) 

whereas, under ex-post pricing we will have: 

v 3 (dj (t + 1)) = 4>j (cj (si (t)) , • • • , 6i (si (t - T))) (32) 

^3,(^ = ^2^ (t). (33) 

Assuming representative agent models, (|30t — d3Tb and (f32T> — (f33T> 
reduce, respectively, to 

v(d(t+l)) = c(d>(d(t),--- ,d(t-T))) (34) 

and 

v(d(t+l))=4>(c(d(t)),--- ,c(d(t-T))). (35) 
Under the persistence model for prediction we have: 

v(d(t + l)) =c(d(t)) (36) 

In the sequel, we will develop a theoretical framework that 
is convenient for analysis of dynamical systems described by 
implicit equations. Such systems arise in many applications which 
incorporate real-time optimization in a feedback loop, several 
instances of which were developed in this section. As we will 
see, this framework is extremely useful for studying the dynamics 
of electricity markets, robustness to disturbances, price stability, 
and volatility under real-time pricing. 

IV. Theoretical Framework 
A. Stability Analysis 

In this section we present several stability criteria based on 
Lyapunov techniques and examine stability properties of the 
clearing price dynamics formulated in Section HTT1 

Theorem 1: Let S be a discrete-time dynamical system de- 
scribed by the state-space equation 

S : x(t + 1) = ip{x{t)) 

(37) 

x e X Q c R+ 

for some function ip : M + M> K + . Then, S is stable if there exists 
a pair of continuously differentiable functions /, g : R + i— >■ M.+ 
satisfying 

g(x(t + l)) = f(x(t)) (38) 
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and 

(i) : 0*=inf{0| f(x) <9\g(x)\, Vxj < 1 (39) 

(ii) : (1l({x I f (x) = g (x)}) = (40) 
and either: 

(iii) : g{x)> 0, Vx, and lim {/ (x) - 5 (x)} < (41) 

x— >oo 

or 

(iii)' : 5 (x) < 0, Vx, and lim {/ (x) - g (x)} > (42) 



Before we proceed with proving Theorem Q] we present the 
following lemma, which will be used several times in this paper. 

Lemma 2: Let X be a subset of R. Suppose that there exists a 
continuously differentiable function / : X H» K, a continuously 
differentiable monotonic function g : X n- M, and a constant 
9 E [0, 00) satisfying 



f(x) <9\g(x)\, WxeX 



Then 



\f(x)-f(y)\<6\g(x)-g(y)\, Vx,yeX 
Furthermore, if (|40t is satisfied, then 

\f(x)-f(y)\ < \g{x)-g{y)\, Vx,yeX, x + y 



(43) 



(44) 



(45) 



Proof: We have 

Vx,yeX, x^y 



domain of ip is K+. Furthermore, as long as / (x (t)) < g (x (t)) , 
the sequence {g (x (t))} decreases strictly. Therefore, (HTt implies 
that 

Vx : 3 M e E, JVeZ + :j(i (t)) < M, Vt > N. (49) 

It follows from ( |49l , monotonicity and continuity of g (•) that the 
sequence {x (t)} is bounded from above too (similar arguments 
prove boundedness of {x (t)} when (|42T > holds). Hence, either 
lim^oo a; (t) = 0, or {x(i)} has a subsequence {x (ti)} which 
converges to a limit x* € R+. In the latter case we have 

lim V (x (t)Him V (x (U)) = lim {/ (x (i,)) - 3 (x (^))} 

= |/(x*)- 5 (x*)| 
If =g{ip(x*)) then c = |/(x*) -g (x*))| = (due to 

H)). If (7(0^(^(0) then 

35, e > 0, s.t. \g{il){x)) -g{x)\ > e, Vx eB(x*,S) 
Define a function r : B (x* , (5) H> R+ according to 

|/(^(x))-/(x)| 
\g^{x))-g{x)\ 

Then it follows from |481 that r (x) < 1 for all x e B(x*,8). 
Furthermore, the function is continuous over the compact set 
B (x* ,5) and achieves its supremum r, where t < 1. Since x (U) 
converges to x* , there exists t € N, such that i(t)e8 (x* , 8) . 
Then 



\f(x)-f(y)\< 



f(r) 



dr 



\g(j)\dT 



where the inequality in d46| i follows from ( |43l and the subsequent 
equality follows from ( |4TT >. Proof of ( |43t is similar, except that 
under the assumptions of the lemma, the non-strict inequality in 
(|46| > can be replaced with a strict inequality. ■ 

We will now present the proof of Theorem Q] 

of Theorem [7} The key idea of the proof is that the function 



V(x (t + l))-rV (x (t)) 

= |/ (x (t + 1)) - / (x (t))\ -r\g(x(t + 1)) - g (x (t))\ 
< 0, Vi > i 

(x) - 3 (y) I (46) Since r < 1. this proves that c = 0. Finally, 

lim /(x(i))= lim <? (x (*)) = <? (x*) 



V(x) = \f(x)-g(x)\ 



(47) 



is strictly monotonically decreasing along the trajectories of d37l ). 
From Lemma [T| we have: 

^(x(t + l))-V(x(t)) 

= |/ (x (t + 1)) - g (x (t + 1))| - |/ (x (t)) - g (x (i))| 
- |/ (x (t + 1)) - / (x (t))| - \g (x (t + l))-g (x (f))| 
< 0. (48) 

Therefore, {V (x (t))} is a strictly decreasing bounded sequence 
and converges to a limit c > 0. We show that c > is not possible. 
Note that the sequence {x (t)} is bounded from below since the 



t— ¥00 t- 
x* = g-H^m f (x (t))) = lim g' 1 o / (x (t)) = lim x (t) 

r— >oo t— j-oo r— ^00 

This completes the proof of convergence for all initial conditions. 
Proof of Lyapunov stability is based on standard arguments in 
proving stability of nonlinear systems (see, e.g., lfL3l ). while using 
the same Lyapunov function defined in ( |47T >. ■ 

Remark 4: The monotonicity conditions in (1411 or (l42l in 
Theorem 03 can be relaxed at the expense of more involved 
technicalities in both the statement of the theorem and its proof. 
As we will see, these assumptions are naturally satisfied in 
applications of interest to this paper. Therefore, we will not bother 
with the technicalities of removing the condition. 

There are situations in which a natural decomposition of 
discrete-time dynamical systems via functions / and g satisfying 
(l38l is readily available. This is often the case for applications that 
involve optimization in a feedback loop, many instances of which 
appeared in section [Ell] For instance, for the price dynamics d25l >. 



we have ip 



and the decomposition is obtained with 
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g = c , and / = v . However, / and g obtained in this way 



may not readily satisfy d39] >. We present the following corollaries. 

Corollary 1: Consider the system (|37| | and suppose that con- 
tinuously differentiable functions /, g : M + h-> R + satisfying 
( 1381 and (|40T> — (f42b are given. Then, the system is stable if there 
exist 9 < 1 and a strictly monotonic, continuously differentiable 
function p : R + H> K satisfying 

p{f{x))f(x)\<6\p{g(x))g{x)\ 

for all x e M.+ . 

Proof: If / and g satisfy ( 1381 ). then so do p o / and p o g 
for any p S C 1 (0, oo). Furthermore, under the assumptions of the 
corollary, pof and pog satisfy d39l— (l42l i. The result then follows 
from Theorem Q] ■ 
Corollary 2: Market Stability I: The system (f25j) is stable 
if there exists a strictly monotonic, continuously differentiable 
function p : M+ n- M satisfying 



for all A 6 



(A)) 



dv- 1 (A) 



9A 



< 



pic 



(A)) 



ac- 1 (A) 



9A 



Similarly, the system (136l l is stable if 

\p (c (x)) c (x) | < 6* |/j (i> (x)) v (x)| 



(50) 



(51) 



for all x S K+. 

Proof: The statements follow from Corollary [T] with f = v^ 1 
and 5 = tr 1 for (l50l >. and f = c and g = v for ( BTl i. and the 
fact that under Assumption I, all of the conditions required in 
Corollary Q] are satisfied. ■ 

Remark 5: By taking p( ) to be the identity function in (T5Qb 
and (IBTl i. we obtain the following sufficient criteria for stability 
of (El or CH: 

|c(x)| < 9\v(x)\ (52) 



cfw- 1 (A) 



9A 



< 



dcr 1 (A) 



o>A 



(53) 



Although these conditions are very simple, they are generally 
harder to satisfy globally for typical supply and demand functions. 

In the economics literature, elasticity is defined as a measure 
of how one variable responds to a change in another variable. In 
particular, price-elasticity of demand is defined as the percentage 
change in the quantity demanded, resulting from one percentage 
change in the price, and is used as a measure of responsiveness, 
or sensitivity of demand to variations in the price. Price-elasticity 
of supply is defined analogously. In this paper, we generalize the 
standard definitions of elasticity as follows. 

Definition 4: Generalized Elasticity: The quantity 



el(X,l) 



A Vdir^X) 



I > 



is the generalized price-elasticity of demand at price A. Similarly, 



el(XJ) 



x Vac- 1 (A) 



(A) 



dX 



I > 



is the generalized price-elasticity of supply at price A. Note that 
these notions depend on the exponent I. For I = 1, we obtain 
the standard notions of elasticity. We define the market's relative 
generalized price-elasticity as the ratio of the generalized price- 
elasticities: 



ii(A,0 = 



(A,0 



(54) 



e p s (A,/) 

The market's maximal relative price-elasticity (MRPE) is defined 

as 

0*(/)= sup |4 (A, 0| ■ (55) 

AGR + 

The notions of generalized demand-elasticity of price and gener- 
alized supply-elasticity of price are defined analogously: 



4(x,o = 



v (x) 



,(x) 



c(x) 



v (x) ' c(x) 
When / = 1, these notions coincide with the Arrow-Pratt 
coefficient of Risk Aversion (RA) J2], |fl6l , and we will adopt 
the same jargon in this paper. The market's relative generalized 
risk aversion factor is defined as: 



RA 

e rel 



(x,l) 



e s p (x,l) c(x) / t>(x) 



e£ (x, /) v (x) \c (x) 



Finally, the market's maximal relative risk-aversion (MRRA) is 
defined as 



v* (0 



sup 



,RA 

~rel 



(x,l) 



(56) 



With a slight abuse of notation, when / = 1, we write e p D (A) 
instead of e v D (A, 1) , and 9* instead of 9* (1), etc. 

The following corollary relates the market's stability to the mar- 
ket's relative price-elasticity (A, /), and relative risk-aversion 

e^(x,l). 

Corollary 3: Market Stability II: The system (J23j is stable 
if the market's MRPE is less than one for some / > 0, that is: 



3/ > 



r (/) = su P | e p el (A, 0| <i 

A 



(57) 



The system (1361 1 is stable if the market's MRRA is less than one 
for some k > 0, that is: 



3Z > : 7]* (0 = sup 



(x,0 



< 1 



(58) 



Proof: The results are obtained by applying Corollary 
criteria (l50l and (|5H . with p (z) = log (z) for I — 1, and 
p{z) = 7 - l+1 for I ^ 1. ■ 

When the cost and value functions are explicitly available, 
conditions (IBTT i or ( 1581 are more convenient to check, whereas, 
when explicit expressions are available for the supply and demand 
functions, it is more convenient to work with (|50t or ((57). 

Example 1: Consider (|25T > with c (x) = x* 3 , and v (x) = 
(x — u) 1 ^ 01 , where a, /3 > 1 and u > is a constant. First, 
consider the u = case. Then, we have 

X(t + l)=P{aX (*))^ 

v (x) = a~ 1 x 5 , v (x) = (1 — a) a~ 2 x 5 

c(x) = /3x^\ c(x) =13(13- l)x ~ 2 
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It can be verified that there cannot exists a constant 9 £ [0, oo) for 
which d52b is satisfied for all x € K+ (equivalently, 9* (0) = oo). 
However, by invoking d58l l with k = 1, we have: 



which satisfy 



< 1 



' \v(x)\ |c(a:)| (a - 1) a- 1 
Hence, the system is stable if 

p < 2 -or 1 



It can be shown that the condition is also necessary and the system 
diverges for j3 > 2 — cT 1 . Moreover, invoking (ISTt with I = 1 
yields exactly the same result, though, this need not be the case 
in general. Consider now the same system with a = f3 = 2 and 
u > 0. Simulations show that the system is not stable in the 
asymptotic sense for u < 1/4. The following table summarizes 
the results of our analysis. 

Table I 





u = 0.25 


u = 0.3 


u = 0.5 


9* (1) = 


2 


2 


2 


9* (1.5) = 


1 


0.872 


0.595 


9* (2) = 


1.299 


0.988 


0.459 



Thus, when u = 1/4, the system is at least marginally stable. 
Furthermore, the above analysis highlights the importance of the 
notion of generalized elasticity introduced earlier (cf. Definition 
IH, as 9* (1) (which is associated with the traditional notion of 
price elasticity) can be greater than one while the system is stable 
and it's stability may be proven using the MRPE for some I > 0. 

The preceding analysis is based on applying the results of 
Theorem[T]and Corollary [3] to systems of the form (f23} (or (|36j), 
which correspond to the persistence prediction model, whether 
it is demand prediction by the ISO in the exante pricing case, 
or price prediction by the consumers in the ex-post pricing case. 
In the next section, we present a theorem that is applicable to 
analysis under the generic prediction models (|26l and ([19V 



< 



g{yk) 



J 71 + 1 



where 



fe=i 



(60) 



(61) 



Then, there exists a constant 70 > 0, which depends only on the 
first n + 1 initial states x (n) , x (0), such that the set 

O = {x e K I 3z e M" : \g (x) -f(x,z)\< 70} (62) 

is invariant under ( |59l , i.e., 

x (T - n) , x (T) 6 fl => x (t) G fid Vi > T 

Furthermore, when (16H holds with strict inequality, the g-scaled 
IAV of x is bounded from above: 

7o 



V s (a0 = 5>(z(t + l))-<?(z(f))| < 



(63) 



fc=i 



Proof: For simplicity and convenience in notation, we prove 
the theorem for the n = 1 case. The proof for the general case is 



entirely analogous. Define the function V : K. 
to 

V(x,z) = \g(x) - f(x,z)\ 

Let 



1 — ^ 



according 
(64) 



70 = V (x (1) , x (0)) + \g (x (1)) -g(x (0))| (65) 

it is sufficient to show 

(66) 



To prove that fl is invariant under 
that 

V(x{T + l),x{T)) < 7o , VT e Z + 



To simplify the notation, define A/t = / {x (t + 1) ,x (t)) — 
/ (a; 0) , x (t - 1)) , and A.g t = g> (i (t + l)j-.g (x (t)) . We have: 

V (x (t + 1) , x (t)) - V (x (t) , a; (t - 1)) 



B. Invariance Analysis 

When functions of the form (|26] i or (fT9l are used for prediction 
of price or demand, the underlying dynamical system is no longer 
a scalar system. An immediate extension of Theorem Q] in its 
full generality to the multidimensional case, while possible, raises 
further complexities in both the proof and the application of the 
theorem. In what follows we take the middle way: we present a 
theorem that exploits the structure of the dynamical system that 
arises from autoregressive prediction models to both make the 
extension possible and to simplify the analysis. 

Theorem 2: Let x : Z + — > R, be a real-valued sequence 
satisfying a state-space equation of the form: 

j(i(t+l)H(x(t),i(i-l),- ,x(t-n)) (59) 

(x (0),...,x(n))EYo Cl" +1 , 

for some continuously differentiable function / : K" +1 h> M, 
and a continuously differentiable monotonic function g : M H> K 



= \g(x (t + l))-f(x (t + l),x(t))\ 

-\g(x(t))-f(x(t),x(t-l))\ 

= \f(x (t),x(t-l))-f(x(t + l),x(t))\ 

-\g(x (t))-g(x (t + l))| 

< \f(x(t),x(t-l))-f(x(t),x(t))\ 

+ \f(x (t),x(t))-f(x (t + l),x(t))\-\Ag t \ 

<9 2 \Ag t _ 1 \ + (9 1 -l)\Ag t \ (67) 

where the first inequality is obtained by applying the triangular 
inequality, and d67l i follows from ( f60b and Lemma|2] By summing 
up both sides of d67l i from t = to t = T we obtain: 

V (x(T + 1) ,x(T)) <V (x(l) ,x(0)) 

T 

+ (0i + 02 - 1) l A ^*l + 62 (l A 5°l " l A ^D ( 68 ) 
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The inequality ( |66l > then follows from ( |68l > and d6Tb . When ( |6TT > 
holds with strict inequality, d63l follows from (l68l l and nonnega- 
tivity of V (x (T + 1) , x (T)) for all Tel. This completes the 
proof. ■ 

It follows from the proof of Theorem |2] that when the initial 
conditions are close to the equilibrium of (|59~l l, it is sufficient to 
satisfy conditions (f60l> — (f6TT> only locally, over a properly defined 
subset of R™ +1 . This is summarized in the following corollary. 

Corollary 4: Let x : Z + -+ R, be a real-valued sequence 
satisfying (|59l l, where / and g are continuously differentiable 
functions. Let 

= {(x, z) G K x K" : \g (x) -f(x,z)\< 70} 

where 70 is given in d64t — (f&Sb . If 

9{Vk) 



oyk 



< 



Vy G fi 



where #fc's satisfy (l6Tb . then Slo is invariant under (|59). Further- 
more, when (loTT i holds with strict inequality, and the initialization 
vector xq — [x (n) ,x (0)] is an element of fio, then (l63l 
holds. 

Theorem|2]and Corollary |4]can be applied to analysis of market 
dynamics under the generic autoregressive prediction models that 
were presented in Section [HI] This includes the generic dynamical 
system models that were developed for price dynamics under 
exante or ex-post pricing (([T8li-(l20ii and (|29j), as well as the 
aggregate demand dynamical systems ( [34-b and (f35T >. The sets 
flo or Hq being invariant implies that the difference between the 
predicted demand and the actual supply (possibly scaled by some 
monotonic function, e.g., log (•)) remains bounded. 

1 ) Analysis of Market Dynamics under Generic Autoregressive 
Prediction Models: In this section we examine the impact of lin- 
ear autoregressive prediction models on market stability. Consider 
the model (l24l >. repeated here for convenience: 

c- 1 (A (t + 1)) = J2l =0 aki- 1 (A (t - k)) 
We apply Theorem |2] (alternatively Corollary |4) with 

g(X) = p(c- 1 (\)) 



(69) 



and 



We examine 



At-n) = P Q 



"'(At-*)) 



-2+1 



\^^k=0 

(|70l i with p (z) = log (z) and p (z) 
I 7^ 1. Conditions (|60l> — (|6TT> then imply that the following 
conditions are sufficient (for some k > 0): 

dw- 1 (A) 

A=A t _ 



Q!fe 



(9A 



E 6 



<9 k e p s (At_ fcj i) 



< 1 



(71) 



(72) 



fe=i 



Conditions (l7Tli-(f72]> are complicated and in general demand 
numerical computation for verification. However, examination of 



(1711 near equilibrium is informative. Suppose that (l24l converges 
to an equilibrium price A. Letting At = Xt-i = ■ ■ ■ = A t _ n = A, 
we observe that the following condition is implied by (FTTb — (T72l>: 



31 > 



k=l 



. (Ik 



: -D 



(A,0|<|e p s (A,OI 



\^k=i 



(73) 



where e p D (A, l) and e p D (A, /) are generalized elasticities as de- 
fined in Definition 0] evaluated at the equilibrium. It can be 
shown that d73l is equivalent to ej? el (A, l) < 1, independently 
of /. Furthermore, for a large class of cost and value functions, 
namely power functions of the form c (x) = x@ and v (x) = x 1 /", 
a,/3 > 1, the equilibrium relative elasticity 9 (A) = e^, (A, l) is 
independent of the autoregressive coefficients ak, k — l,..,n. 
Thus, if the closed-loop market is unstable under the persistent 
prediction model (ai = 1, ak = 0, k ^ 1), then global stability 
cannot be verified for any linear auto-regressive model of the form 
d24l) using (TTTb— (T72b. Although this analysis is based on sufficient 
criteria, it suggests that it may be difficult to globally stabilize 
these systems via linear autoregressive prediction. Indeed, exten- 
sive simulations show that such models will not globally stabilize 
an unstable market, unless the MRPE is very close to one. For 
values of 8* > 1.05 global stabilization could not be achieved 
in our simulations. Local stabilization is, however, possible for 
moderate values of 6*, namely, 9* ^ 3. 

2) Analysis of Dynamics of Markets with Exogenous Inputs: 

Theorem 3: Let x : Z+ — > R and u : Z+ -+ M be real-valued 
sequences which satisfy a state-space equation of the form: 



g(x(t + l))=f(x{t),u{t)), u(t)eU 

x (0^o C M 

for some continuously differentiable function / : M 2 
a continuously differentiable monotonic function g : 
satisfying 



8 it \ 



< 1, Vx e R, u g U 



and 



9 ft \ 



< 



g{x) 



, Vx G R, u G U, 



(74) 

1 and 

-> R 

(75) 
(76) 



(70) where 



U = {11 G M : |u| < k} 

and k G (0, 00) , and 8 G [0, 1). Define 

1 + 9 



C« (0) = « 



1 



(77) 



Then, the set 

fi(fi) = {a::||/(a,i/)-ff(a:)|-|i/|| <C«(fl), G C/} (78) 

is invariant under (l74i i. Furthermore, the 5-scaled IMV of x is 
bounded from above: 



i — »-oo J f 1 



2k 



(79) 
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Proof: Define 
V(x) = su V { \f{x,v)-g{x)\-\v\ } - £ K (6) . 

It is sufficient to show that there exists r > 0, such that: 

V(x(t + 1)) -rV{x{t)) < 0, Vt e Z+. 

To simplify the notation, define Aft = f (x (t + 1) , u (t + 1)) — 

/ (x (t) , u (t)) , and A. 9t = g (a: (t + 1)) - <? (x (<)) . Then 

V(a; (t + l))-rV (a:(i)) 

sup{||/(s (t+l),i/)- ff (a: (t + 1))| - | } 

{\\f{x{t),v)-g{x{t))\-\v\\) 

+ Ck(8) (r-l) 



t sup 

ueu 



<SU P f |/(x(t + l), !/)-/(!(*), «(t))|-H } 

-r|A St |+r/s + C.W(r-l) (80) 



< sup|/(a:(f + (a; (t),i/)| 

!/££/ 



+ SU P { |/ (*(*),!/)-/ (3 (t),tt(t))|-|l/| } 

-T\Ag t \+TK + ( K (6)(T-l) (81) 
< (0 - r) |Aflt| + (1 + r) K + Ck (0) (r - 1) (82) 

where (l80T > follows from the choice of v = u (t) and |it (i) < k, 
(IHTT i follows from the triangular inequality, and (l82t follows from 
(I75t -(l761l and Lemma [2] The desired result follows from the fact 
that the right-hand side of (l82l will be non-positive for t = 0, 
and Q (0) defined in (|77). To prove <[79), let r = 1 in (f82l to 
obtain 

V(x(t+ 1)) - V (a; (*)) < {9 - 1) |Ag t | + 2k (83) 
Summing both sides of (l83l over all t = 0, 1, .., T results in: 

T 

V (x (T + 1)) < V (a (0)) + (0 - 1) |Aff t | + 2Tk (84) 

t=i 

It follows from d84) and non-negativity of V (x (T + 1)) + ( K (9) 
that 

T 

(1 - 9) |Aff t | < 2Tk + y (x (0)) + C« (9) • (85) 

t=i 

The desired result (l79l then follows immediately from (l84l by 
dividing by T and taking the limit as T — > oo. ■ 
The following corollary is a local variant of Theorem [3] and 
is useful for scenarios in which, there exists no positive number 
9 < 1 such that d76b is satisfied for all whereas it might 

be possible to satisfy the inequality locally over a subset that 
contains ft (9). 



Corollary 5: Let x : Z+ — > R and u : Z + — ► M. be real-valued 
sequences satisfying d74b . For 6* < 1, define: 



= mi{9 



0_ 

dx 



f{x,u) 



< 



d_ 

dx 



Vx G O (6/) , «£(/ 



where fi (0) is given in ( f78l . Then f2 f^*J i s i nvar i ant under (|74] > 
if 

6»* < 6». 

Furthermore, ((79) holds with 9 = 9*. 

Consider equation d74| i or (|59) . When the functions g and / are 
p-scaled supply and demand functions, the minimal 9 satisfying 
d76T i or (f60b will be the MRPE associated with the market models 
(T74b or (|59] l. When g and / are p-scaled marginal value and 
marginal cost functions respectively, the minimal 9 satisfying the 
inequalities will be the MRRA associated with d74l ) or d59l ). In the 
remainder of this section, we consider applications of Theorem 
[3] to the two time-varying models of consumer behavior (fl4l and 

aj Multiplicative Perturbation: Consider the multiplicative 
perturbation model (fT2t . Under this model, the market dynamics 
is given by 

c- 1 (A (t + 1)) = (l + Is (t)\ v- 1 (A (t)) , 5 (t) e [-«, «] 

(86) 

where the 1/2 factor in front of 8 (t) is simply a scaling factor. 
We invoke Theorem [3] with 



5 (A) = log (c- 1 (A)) 



(87) 



and 



f(X,S)=log(l + S/2)v- 1 (X) 

= log (1 + 5/2) + log Or 1 (A)). 

It can be verified that d75l l and d76l l are satisfied as long as n < 1 
and 9* < 1, where 6 1 * is the MRPE defined in 05b . Furthermore, 
( K (9* ) is the upperbound on the size of the invariant set, where 
Ck (•) is defined in ( f77l i. In particular as 9* — > 1, small perturba- 
tions may induce extremely large fluctuations as measured by log- 
scaled IMV of supply. The theoretical upperbound is 1/ (1 — 9* ) . 
When Corollary [5] is applicable, the size of the invariant set can 
be characterized by Ck(#*), where 9* is the market's local relative 
price-elasticity. Furthermore, volatility can be characterized by 9* 
as well. 

b) Additive Perturbation: Under the additive perturbation 
model (fl4l . the market dynamics can be written as 

(T 1 (A (t + 1)) = M + ]rU (t) + V- 1 (A (£)) , U{t) G [-K, K} 

(88) 

where m > 1 is a shifting factor, and k < u , so that the demand 
is always at least uq/2. We invoke Theorem [3] with d87l i and 

/ (A, u) = log (u Q + + v^ 1 (A) 
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Then, under the above assumptions, d75l l is satisfied. In a similar 
fashion to previous analyzes, (l76l > can be related to the market's 
relative price-elasticity. In this case, the price-elasticity of demand 
turns out to be: 

df{X,u) X flu" 1 (A) 



£d(A) 



dX 



ox 



u a + u/2 + v- 1 (A) 

The larger the minimum of the inelastic component (i.e., no — 
k/2), the smaller the price-elasticity of the overall demand will 
be. Under the assumptions made above, there is always a nonzero 
minimal demand u m i n (i) — Uq/2. Therefore, it is sufficient to 
verify ( TT5T > over A > c (u /2) instead of all A > 0. In conclusion, 
(|76l i reduces to: 



dv- 1 (A) /OX 



uq/2 + v- 1 (A) 



< 



dc- 1 (A) /8X 



-i(A) 



VA > c (uq/2) (89) 



Let 9* be the minimal 9 satisfying ( |89l ). Similar to the case with 
multiplicative uncertainty, in this case too, the upperbound on the 
size of the invariant set is given by where Cre (') i s given 

in d77b . Moreover, the log-scaled IMV of supply is upperbounded 

byuo/(l-fl*). 

The analysis reconfirms the intuition that participation of a 
small portion of the population in real-time pricing will not have 
a severe destabilizing effect on the system, as satisfying d89b for 
larger values of uq is easier. System stability concerns should 
arise when a large portion of the population is exposed to real- 
time pricing. 

Remark 6: It can be proven that when u (t) is a periodic 
function with period T and d76l l is satisfied, then all solutions 
of d88l converge to a periodic trajectory with period T, 



C. Volatility 

The following corollaries follow from Theorems [2] and [3] and 
explicitly relate the market's MRPE and MRRA to volatility. 

Corollary 6: Volatility I: Let 9* < 1 and r)* < 1 be the MRPE 
and MRRA associated with the market model ( 1741 ). Then, there 
exists a constant C, depending on the size of the disturbances 
only, such that the log-scaled IMV of supply is upperbounded by 
C/ (1-0*), i.e., 

1 T 

lim llogtc" 1 (A(t+ 1))) 



And the log-scaled 
C/(l-rf),i.e., 



logOT 1 ^))) 



< 



c 



IMV of price is 



1 - 9* 
(90) 

upperbounded by 



lim -Y 

t=l 



log(A(t+l))-log(A(i))|< 



C 



1 



T) 



(91) 



Corollary 7: Volatility II: Let 9* < 1 and rf < 1 be the 
MRPE and MRRA associated with the market model < f59b with 
linear autoregressive prediction. Then, there exists a constant 
C such that the log-scaled IAV of supply is upperbounded by 
C/ (1 - 6*) , i.e., 

c 



E 

t=l 



log 



- 1 (a (t + 1))) - log (cr 1 (a (t)))| 



< 



i 



(92) 



And the log-scaled IAV of price is upperbounded by C/ (1 — rf) , 
i.e., 

oo „ 

^|log(A(t + l))-log(A(t))|< r 



(93) 



Remark 7: Generalized versions of the above corollaries can 
be formulated based on 9* (I) and -q* (I) , in which case the 
scalings of the signals need to be defined accordingly: letting 
pi (x) = x~ l+1 for 1^1, the p;-scaled IMV of supply and price 
will be upperbounded by C/ (1 - 9* (I)) and Cj (1 - r/* (I)) 
respectively. Furthermore, when the prices remain bounded within 
an invariant set, e.g., when the conditions of Corollary 2] or 
Corollary [5] hold, one can replace 9* (I) and rj* (I) with local 
relative elasticity ratios 9* (I) and rj* (I) . 

D. Robustness and Incremental L2-Gain 

The p-scaled incremental L2-gain of a discrete-time dynamical 
system with input signal u : Z — > K and output signal h : Z — > K 
is defined to be the minimal 7 > such that the inequality 



7||p(u)-p(«)|| 2 -||p(fe)-p(&)|| 3 >0 



(94) 



is satisfied for all input/output pairs (u, h) and (u, h) such that 

p(u) - p{u) € £2- 

For systems with larger gains, it is generally expected that 
relatively small deviations from a nominal input u would stir 
relatively larger deviations from the nominal output signal h. 
This gain can be used as a metric for assessing the robust- 
ness/sensitivity of the system to arbitrary external disturbances. 
It can be proven that for the market model (l86t (more generally, 
the market model (l74l i with multiplicative uncertainty), the log- 
scaled incremental L2-gain from the perturbation 8 (•) to the 
demand is upperbounded by 9* / (1 — 9*) . The gain from i5 (•) to 
the supply is upperbounded by 1/ (1 — 9*) . These results-stated 
here without proof-quantify the dependence of the closed-loop 
system's robustness, as measured by the incremental L2-gain, on 
the markets maximal relative price-elasticity. 

V. Discussion 

Cho and Meyn [7] have investigated the problem of volatility 
of power markets in a dynamic general equilibrium framework. 
Their model can be viewed as a full-information model in which 
the system operator has full information about the cost and value 
functions of the producers and consumers. Market clearing is 
instantaneous and supply and demand are matched with no time 
lag. The producer's problem is, however, subject to supply friction 
or a ramp constraint, i.e., a finite bound on the rate of change of 
supply capacity. It is concluded that efficient equilibria are volatile 
and volatility is attributed to the supply friction. In the formulation 
of the consumer's problem is not subject to ramp constraints. 
In our formulation, neither the consumer's problem nor the 
producer's is explicitly subject to ramp constraints, yet other 
factors are shown to contribute to volatility, namely, information 
asymmetry and high price elasticity of demand. Interestingly, 
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if we included ramp constraints in the consumer's problem it 
would have a stabilizing effect, as it would limit the consumer's 
responsiveness to price signals and reduce her elasticity. This 
effect is implicitly and qualitatively captured in our framework 
through the introduction of an inelastic component in the demand, 
which certainly limits the rate of change in the demand, and 
was shown to have a stabilizing effect. However, uncertainty 
in the supply side, either in the available capacity or in the 
cost, works in the reverse direction: when supply is sufficiently 
volatile, a trade-off might exist and responsiveness and increased 
elasticity of demand might be desirable, though this needs to be 
quantified rigorously. The models developed in the paper do not 
include uncertainty in generation, and this would be an interesting 
direction for future research. 

The above discussion leads to another interesting question: 
"quantifying the value of information in closed-loop electricity 
markets". Given the heterogeneous nature of consumers and 
time-varying uncertainty in their preferences, needs, and valua- 
tions for electricity, learning their value functions and predicting 
their response to a price signal in real-time appears to be a 
difficult problem. Suppose that the consumers provide a real- 
time estimate of their inelastic and elastic consumption to the 
ISO. How valuable will this real-time information be and what 
would be its impact on volatility and reliability of the system? 
Given the potentially significant costs and barriers associated with 
obtaining such information in real-time, quantifying the value of 
information in this context seems an extremely important and 
timely question with potentially significant impact the architecture 
of future power grids. 

VI. Numerical Simulations 

In this section we present the results of some numerical 
simulation. For the purpose of simulations, we use the following 
demand model: 

D (t) = fxidt (t) +fi 2 (l + 6 2 (t)) v- 1 (A (<)) (95) 

where d\ (t) is the exogenous, inelastic demand: 

di (t) = ciq + ai sin (t) + a 2 sin (2t) + S\ (t) 

and Si (t) - 7V(0,0.1 2 ) and S 2 (t) ~ TV (0, 0.01 2 ) are random 
disturbances. The parameters /ii and /i 2 are adjusted, on a case-to- 
case basis, such that the average demand under real-time pricing 
(i.e., when [i 2 > 0, fi% < 1) remains nearly equal to the average 
demand in the open loop market (p, 2 = 0, /ii = 1), that is: 

This normalization, takes out the effect of higher or lower average 
demand on price and allows for a fair comparison of volatility 
of prices in open-loop and closed-loop markets. The following 
parameters are chosen for all simulations in this section: 

a = 4 GW, ax = 1 GW, a 2 = 1 GW 

This puts the peak of the inelastic demand at 6 GW and the valley 
at 2 GW, modulo the random disturbance 8\ (t) . All simulations 



are for a 24 hour period and prices are updated every 5 minutes. 
The average demand in all simulations is approximately 4 GW 
per five minutes for both open-loop and closed-loop markets. The 
metric for comparison in these simulation is the Relative Volatility 
Ratio (RVR), defined as the ratio of the log-scaled IAV of the 
closed-loop market to the log-scaled IAV of the open-loop market. 
The results of the first simulation are summarized in Figure 
The prices are extremely volatile under real-time pricing (RVR 
= 51.12) and the system is practically unstable. 
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Fig. 3. Simulation of a market with quadratic cost function c(x) = x 3 , value 
function v (x) = log (x) , and demand function D (t) given in (95) with /ii = 
0.075, fi 2 = 2. 

The results of the second simulation are summarized in Figure 
[4] Based on the chosen parameters, this market is less volatile than 
the one in the first simulation, yet, volatility of demand increases 
under real-time pricing (RVR=2.33). Since in this simulation the 
cost is quadratic, the price (not shown) has a very similar pattern. 
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Fig. 4. Simulation of a market with quadratic cost function c(x) = 3a; 2 , 
value function v (x) = y/x, and demand function D (t) in (95) with /ii = 0.7, 
fi 2 = 3 X 10 3 . 

The third simulation is summarized in Figure |4] For each 
value of /ii G [0, 1] (with 0.05 increments), the expected RVR 
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was calculated by taking the average RVR of 50 randomized 
simulations. The random parameters are 61 (t) , 62 (t) , and the 
initial conditions. The experiment was repeated for four different 
value functions: v (x) = x x l a , a — 4,4.5,5,5.5. It is observed 
that volatility increases with decreasing a or /ii, both of which 
increase the price-elasticity of demand. 



Relative Volatility versus Penetration of Real-Time Pricing 
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Fig. 5. Simulation of a market with quadratic cost function c(x) = 3a; 2 , value 
function v(x) = x 1 /" , and demand function D(t) given in __} with /^i £ [0, 1], 
and /i2 adjusted accordingly to keep the total demand constant. Decreasing a or fi\ 
increase the price-elasticity of the overall demand and hence, increase volatility. 



VII. Conclusions and Future Work 

We investigated the effects of real-time pricing on the stability 
and volatility of electricity markets, and showed that exposing the 
retail consumers to the real-time wholesale market prices creates a 
closed-loop feedback system which could be very volatile or even 
unstable. When the system is stable, an upper bound on volatility 
and robustness to external disturbances can be characterized in 
terms of the market's relative price-elasticity, defined as the 
ratio of generalized price-elasticity of consumers to that of the 
producers. As this ratio increases, the system may become more 
volatile, eventually becoming unstable when the ratio exceeds 
one. As the penetration of new demand response technologies and 
distributed storage within the power grid increases, so does the 
price-elasticity of demand, and this is likely to increase volatility 
and possibly destabilize the system under current market and 
system operation practices. While the system can be stabilized 
and volatility can be reduced in many different ways, e.g., via 
static or dynamic controllers regulating the interaction of whole- 
sale markets and retail consumers, different pricing mechanisms 
pose different limitations on competing factors of interest. In 
light of this, systematic analysis of the implications of different 
pricing mechanisms, and quantifying the value of information 
and characterization of the fundamental trade-offs between price 
volatility and economic efficiency, as well as system reliability 
and environmental efficiency are important directions of future 
research. In summary, more sophisticated models of demand, a 



deeper understanding of consumer behavior in response to real- 
time prices, and a thorough understanding of the implications of 
different market mechanisms and system architectures are needed 
before real-time pricing can be implemented in large-scale. 
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